Finite Temperature Density Instability at High Landau Level Occupancy 
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We study here the onset of charge density wave instabilities in quantum Hall systems at finite 
temperature for Landau level filling v > 4. Specific emphasis is placed on the role of disorder as well 
as an in-plane magnetic held. Beyond some critical value, disorder is observed to suppress the charge 
density wave melting temperature to zero. In addition, we find that a transition from perpendicular to 
parallel stripes (relative to the in-plane magnetic held) exists when the electron gas thickness exceeds 
~ 60A. The perpendicular alignment of the stripes is in agreement with the experimental finding that 
the easy conduction direction is perpendicular to the in-plane held. 



Recent experiments by Lilly and co-workers [lj] as well 
as by Du and co-workers |^| have shown surprisingly that 
when a 2D electron gas is subjected to perpendicular 
magnetic fields of intermediate strength so that multiple 
Landau levels (N > 2) are occupied, the resultant trans- 
port properties are strongly anisotropic. Specifically, the 
longitudinal resistivity was found to depend on the di- 
rection of the current at sufficiently low temperature for 
filling factors v > 4. While the anisotropy is largest at 
v = 9/2 and decreases subsequently as v is increased, 
the signature anisotropy is observed to persist for filling- 
factors up to v — 23/2. This behaviour is in contrast 
to what is traditionally observed at high magnetic fields 
such that only the N = or N = 1 Landau levels are oc- 
cupied |p|-p|. In addition, subsequent experiments have 
shown [|6|,[7| that the axis along which the anisotropy is 
oriented can be changed by a sufficiently strong in-plane 
magnetic field. Further, at v — 5/2, an anisotropic state 
is observed only in the presence of a parallel field. This 
indicates that in the presence of an in-plane field, the 
anisotropic electronic state is lower in energy than the 
corresponding isotropic state at v = 5/2. 

Anisotropic transport is inconsistent with a uniform 
electron state. As a consequence, charge density wave 
(CDW) formation has risen to the fore as the leading 
explanation of these experiments [jl]Q] . In principle, the 
ground state of a 2D electron gas is expected to be non- 
uniform when the magnetic field (£?) is sufficiently low so 
that no fractional quantum Hall states are accessible but 
not so low that disorder destroys the gap established by 
the cyclotron frequency, lo c . In this intermediate range 
of magnetic fields, the ground state is governed by the 
competition between the short-range exchange and long- 
range direct Coulomb interactions. Because these in- 
teractions are of opposite sign, the electron gas breaks 
into domains. Several years ago, Fogler, Koulakov, and 
Shklovskii @(FKS) and Moessner and Chalker § (MC) 
showed that at the Hartree-Fock level at T = 0, the filling 
in the highest Landau level oscillates between zero and 
unity with a well-defined period. Periodic stripe or bub- 



ble phases were found to be stable as the filling in the top 
Landau level was varied. Recent numerical simulations 
on small numbers of particles have also confirmed 
charge density wave formation for N > 2. 

In this work, we focus on the role disorder and a fi- 
nite in-plane magnetic field play on the CDW instabil- 
ity at finite temperature for N > 2. First, we include 
disorder in an iterative screened Hartree-Fock approach 
that incorporates Landau level mixing explicitly. Second, 
we compare our results with those of previous effective 
theories for the top Landau level. We find that 

the results from the effective theory for the width of the 
stripes and the melting temperature are in good agree- 
ment with the approach used here. We then adopt the 
effective interaction approach to study the role of an in- 
plane magnetic field. The key assumption is that it is 
the orbital effect of the in-plane field rather than the 
Zeeman energy that affects the CDW instability. In ad- 
dition, we find that a transition from perpendicular to 
parallel stripes (relative to the in-plane magnetic field) 
exists when the electron gas thickness exceeds ~ 60A. 

The starting point for our analysis is the Hartree-Fock 
(HF) equation, (H + T^ F )ip 

where Hq is the non-interacting Hamiltonian for a sys- 
tem of electrons in a transverse magnetic field, is 
the Hartree-Fock self-energy, %b naxo (r) are the single par- 
ticle wavefunctions and E niJXo are the single particle en- 
ergies. This matrix equation leads naturally to Landau- 
level mixing when the self energy, E HF , is non-diagonal 
as is the case for an inhomogeneous system. To solve this 
equation, we expand the single particle wavefunctions 

(r) = ^ C nn i (a, xo)cf) n > Xo (1) 

n' 

in terms of the Landau basis, 
?W ( r ) = ex P(- ix oy/£ 2 )® n (x - x )/y/L^, where $„ is 
a harmonic oscillator state centered at xq = 2ir£ 2 p/L y 
with p — 0,±1, • • • and L y is the width of the sample 
in the y-direction. Due to computational constraints, we 
limit our analysis to the configurations with translational 
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symmetry in one direction. Therefore we are concerned 
only with stripe-like instabilities. In the Landau basis, 
the HF equations 



^mn'Crm'f^ ^o) ^iutxq Cnm (^") ^o) 



(2) 



reduce to a self-consistent determination of the expansion 
coefficients C nn > and E naxg . The matrix elements of the 
Hamiltonian involve a free term, (Ho) nn i — 5 nn i(Tiu> c (n-\- 
1/2) — agfi^B /2), where \xb is the Bohr magneton and 
g the gyromagnctic ratio as well as contributions from 
direct and exchange self-energy interactions: 



mayo ) (4 > nx ' t Pmayo \^\^n' x ^Pmay ) 

mayo 

^nn'( CTX o) = ~ ^ f{E mcr y )((/) nXa 1p m<7 y \u\ljj m <Ty <l>n'xo) 
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. In the di- 

= 2^e 2 / K | g |, 



with f{E nayo ) = [exp((E ntTyo - (jl)/T) + 1] 
rect term, the Coulomb interaction, v(q) 
is screened only by the semiconductor background with 
static dielectric constant k, while the exchange term, 
u(q) — 2ne 2 / ne(q)\q\, is screened by the electron liquid 
with dielectric constant, e(q) = 1 — (2ne 2 / q)x(q) [H - 

At the level of the random-phase approximation 
(RPA), the exact expressions for this polarization effect 
can be calculated only for a homogeneous system. To ex- 
tend these results to the inhomogeneous case, we perform 
an average of x over the period of the CDW. This proce- 
dure should be accurate as long as x is slowly varying in 
the range of filling factors spanned by the period of the 
CDW. Consequently, following the result of Manolescu 
and Gerhardts Jl2[], we write the inter-Landau suscepti- 
bility as 
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where R is the number of guiding centers within a period, 



F nn - (Z) = -f ) z (n-n')/2 e -z/2 L n-n' {z) (4) 



and L™(z) is an associated Laguerre polynomial. 

In the presence of disorder, the energy levels acquire 
finite width. Experimentally, the broadening parameter, 
r, has been found to scale as T = Tq^/B. Phenomenolog- 
ically, this behavior can be incorporated into our scheme 
by using the spectral function in the form [O 



in 

P*(n,yo,uj) = —=-exp 
V 71 " i 



(huj — E nayo ) 21 
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(5) 



As a consequence the Fermi functions f(E nayo ) in the 
above equations need to be replaced by the filling factors 

duf(huj)p a (n,y ;u;). (6) 



We solve Equations in real space starting from 

a non-interacting initial guess for the energy levels. To 
evaluate the inhomogeneous screening, we chose a grid 
composed of 16 to 32 points per period of the CDW. We 
included as many as 10 Landau levels. The parameters 
used in the computation (which are specialized to GaAs) 
are as follows: m* — 0.067m e , g — 0.4, k = 12.7, and 
n = 2.4 x 10 1:l cto~ 2 . With these parameters, Huj c = 20-B, 
for the cyclotron gap, e 2 / n£ = for the strength of 

Coulomb interaction, and 0.133-B for bare spin splitting. 
Here, the energies are measured in Kelvin and magnetic 
field, B, in Tesla. All of the computations were performed 
for a fixed number of electrons, with the chemical poten- 
tial, fx, being determined self-consistently. 

In order to obtain the CDW versus homogeneous phase 
diagram, we must first determine the optimal period of 
the CDW. The optimal period yields the highest melting 
temperature for a given average filling. We find that this 
wavelength coincides with the period of the CDW with 
the lowest free energy. The results for the optimal period, 
Ai, at half filling in terms of magnetic length I are shown 
below: 



V 


4.5 


5.5 


6.5 


7.5 


Ai 


5.84 


5.90 


6.93 


6.99 


A 2 


6.13 


6.19 


7.18 


7.22 


A 3 


6.04 


6.04 


7.14 


7.14 



The error in the determined values is on the order of 
1%. From the theory of FKS, it is expected that the 
period at half- filling, A3 , of the CDW should only depend 
on N. Indeed, we observe minor variation for different 
spin polarizations, and the periods we obtain are quite 
close to the FKS value, A 3 = 2.7y/2N+l. Within the 
same Landau level, we find however, a reduction in A 
away from half-filling. For example, at filling fractions 
0.2 and 0.8, the reduction is on the order of 4 — 5%. The 
values of A2 were obtained from an effective theory for 
the top Landau level 0. These values are also in close 
agreement with the theoretical estimate and the results 
from the numerical scheme outlined above. 

Shown in Fig. (^) are the melting temperatures at half 
filling as a function of magnetic field. The triangles repre- 
sent the results for the screened Hartree-Fock approach. 
As mean-field estimates these values will certainly be re- 
duced by quantum fluctuations Jl4] | and hence set an 
upper bound for the temperature at which a charge in- 
stability occurs. These values do, however, conform to 
the experimental trend that the temperature at which 
the anisotropy is observed decreases as v increases. Also 
shown in Fig. (^) are the melting temperatures computed 
using the effective theory B for the top Landau level. 
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The circles represent melting temperatures at half-filling 
and the asterisks to the corresponding quantity for an 
electron gas with root-mean-square (rms) thickness of 
100 A. As is evident, all three estimates for the CDW 
melting temperatures are in close agreement suggesting 
that the approximations incurred in the effective theories 
are consistent. The inset in Fig. ([[]) shows how the dis- 
order as measured by r attenuates CDW formation at 
v = 4.5. As is evident, a critical value of the disorder, 
rg s» 8.3, exists above which the CDW phases desist. 
Experimentally, the disorder parameter Tq can be esti- 
mated, for example, from the lowest magnetic field for 
which de Haas- Van Alphen oscillations are observed. In 
very clean samples jj], the oscillations disappear below 
50 mT. Assuming that the oscillations disappear when 
2r ~ fuj c , we estimate that T ~ 2.2[K /T 1 / 2 ]. We find, 
that for the range of fillings studied (y < 15/2) the crit- 
ical value To > 2.2. 
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FIG. 1. Melting temperatures for CDW formation at 
half-filling as a function of the perpendicular magnetic field 
for three cases: 1) A represent the melting temperatures from 
the screened Hartree-Fock approach, 2) o were calculated from 
Eq. ([?]) for zero thickness of the 2D layer, and 3) * the same 
method but with an rms thickness of 100A. The inset shows 
the melting temperature as a function of disorder for v = 4.5. 



To understand the role of an in-plane field, we adopt 
the simpler approach in which the interactions in the top 
level are treated within an effective theory fslJslJllll . This 
simplification is warranted as we have seen that the re- 
sults of the effective theories are in agreement with the 
more elaborate scheme described above. In the spirit 
of the work of Fukuyama, Platzman and Anderson [fl5| , 
we work within an effective Ginsburg-Landau theory and 
write the melting temperature for a CDW in a 2D elec- 
tron gas in a perpendicular magnetic field as 

T c = v(l-v)\U(Q )l (7) 

with v is the fractional occupancy in the top Landau 



level and U(Qo) is the effective Hartree-Fock potential 
evaluated at the momentum Q where the potential is 
minimized. 

The primary hurdle in computing T c is the effective 
Hartree-Fock potential U(Qq). As the spin degrees of 
freedom do not couple to the stripe orientation, the spin- 
splitting induced by the in-plane field cannot account for 
the experimental observations. We focus then on the 
orbital effect of the in-plane field. To this end, we con- 
sider our quantum Hall system to have a finite thickness 
in the direction normal to the 2D plane. Due to this 
finite thickness, the one particle wave functions couple 
to the in-plane component of the magnetic field Bn. In 
what follows, we will neglect the contribution due to the 
coupling of £?|| to the electron spin. Let us consider a 
system of electrons moving in the x — y plane in the 
presence of a tilted magnetic field B = (Bt&n9, 0,B). 
We assume that the confining potential in the z direc- 



tion is harmonic V(z) 



*fl z /2, with a character- 



istic frequency determined from the rms thickness in 
the z-direction, = h/2mL 2 . While the actual con- 
fining potential differs from the parabolic form chosen 
here, the main features of the finite thickness should 
be correctly captured by appropriately choosing the rms 
width. By solving the Schrodinger equation in the gauge 
A = (0, — zBtm\9 + xB,0), we obtain the one-particle 
wave functions for this system: 

0nx o (r) = —7=e $~$q + ((x - x )sm8 + zcosO) 



x ((x - Xq) cos( 



(8) 



where $^ is the harmonic oscillator wave function cor- 
responding to the frequency to and 8 is given by tan 9 — 
t&nOtu 2 / {tjj\ — uj 2 ). The frequencies u>± 



4 - \^ 



1 



+ tt 2 u 2 tan 2 9 



depend on the strength and orientation of the magnetic 
field and on the confining potential. In the absence of 
Bu, uj+ reduces to and lj_ to the cyclotron frequency 

LU C . 

Using the one-particle states Eq. (B), we are able now 



to compute the matrix elements 



3 iQ-r| 



needed 



for the determination of the effective potential U(Qo). 
Explicitly, we write the effective potential as a sum of a 
Hartree term 



ne{q Xl q y )ll 
and a Fock term 



where F% n > (q) is the analog of Eq. (0) for the case of a 
tilted field and is given by 
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with 

7 2 = (g, cos + q x sin tf) 2 ^ + g 2 sin 2 Utjl\ 

a 2 = (q x cos 6 - q z sin 6) 2 l 2 _ + q 2 y cos 2 6ll/l 2 _. (9) 

The dielectric constant e(q x ,q y ) is calculated in the RPA 
approximation (see Eq. (0)). In this case, we used the 
non-interacting energies E® a — huj c (n + i) — ^gfisB 
instead of the self-consistent energies E naXQ and the q 
dependence enters through the angle dependent function 

Using this method, we minimized the effective poten- 
tial for two orientations of the stripes: perpendicular 
and parallel to the in-plane magnetic field. In general, 
the minimum value of the effective potential is displaced 
both vertically in energy and horizontally in wavevector. 
With extreme consistency, we observed that the wavevec- 
tor for stripes oriented parallel to the in-plane field always 
exceeded the wavevector for those perpendicular to the 
stripes. This difference was not more than 4%. Our re- 
sults for the change in the melting temperature (T c ) as 
a function of an in-plane field for an rms thickness of 
fOOAare shown in Fig. (0). The solid curves correspond 
to a parallel orientation of the stripes (with respect to the 
in-plane field) while the dashed curves denote the per- 
pendicular orientation of the stripes. As is evident, the 
melting temperature for stripes oriented parallel to the 
field is consistently higher than T c for the perpendicular 
orientation as v is incremented. However, as shown in the 
inset, this trend reverses as the rms thickness is reduced 
below 65.3Afor v = 6.5. Below this thickness, the differ- 
ence in melting temperature between perpendicular and 
parallel stripes is 0(lmK) for a magnetic field tilted at 
45°. This temperature difference directly translates into 
an energy difference per electron between the two phases 
at low temperatures || . To determine which phase is pre- 
ferred, we must compare the total energy difference with 
the temperature. As the temperature is lowered below 
T c , the correlation length of the ordered charge density 
wave domains grows until perfect orientational order is 
reached at the Kosterlitz-Thouless temperature. For suf- 
ficiently large domains, even miniscule energy differences 
per electron can lead to orientational order and hence 
anisotropic transport. Hence, we predict that for thin 
electron gases (<60A), the perpendicular orientation is 
preferred, whereas for thicker electron gases (>60A) par- 
allel stripes dominate. Interestingly enough, Eisenstein 
p6[ places the electron gas thickness at 58A and experi- 
ments identify the direction of the in-plane field as be- 
ing the high resistive direction. This experimental result 
is consistent with our finding that perpendicular stripes 
are energetically favoured for thin electron gases. 
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FIG. 2. Change in melting temperature of the CDW stripe 
phases as a result of a parallel magnetic field for an rms thick- 
ness of fOOA. The solid curves correspond to stripes oriented 
parallel to the field and the dashed lines to a perpendicular 
alignment with respect to the field. Four filling factors are 
shown: □, o, O, and A correspond to v =4.5, 5.5, 6.5, and 
7.5, respectively. The inset shows the difference in the melt- 
ing temperature between parallel and perpendicular stripes 
as a function of the thickness of the electron gas in A. 
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